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1.  Background 


Nuclear  batteries  ean  have  a  small  but  signilieant  impaet  as  a  power  souree  for  long-lived 
sensors  {1-3).  For  example,  eonsider  embedded  sensors  for  bridges  and  foundations  of  large 
buildings  where  these  sensors  ean  provide  information  about  the  struetural  integrity  {4,5). 
Unattended  sensors  for  oeeanie  temperature  measurement  would  need  no  maintenanee  during  an 
increased  operational  lifetime.  Unique  opportunities  are  waiting  for  sensor  applieations  that 
would  benefit  from  keep-alive  power  lasting  30  years. 

Chemieal  power  sourees  ean  store  the  total  energy  required,  but  eleetrieal  leakage  and 
deterioration  stops  them  from  sustaining  the  requisite  power  over  long  lifetimes.  Nuelear 
batteries  offer  a  unique  souree  to  fill  the  niehe.  Even  though  the  energy  storage  density  of 
nuelear  isotopes  is  six  orders  of  magnitude  greater  than  that  possible  in  chemieal  bonds,  the 
stigma  assoeiated  with  nuelear  materials  will  require  eareful  design  and  thoughtful  applieation. 
Sensors  embedded  in  structures,  providing  information  useful  to  monitor  struetural  integrity  may 
be  just  the  applieation  to  overeome.  The  negatives  often  assoeiated  with  nuelear  applieations. 

MCNPX  is  a  general-purpose  Monte  Carlo  N-Partiele  eode  whieh  supports  34  partiele  types, 
eontains  cross-seetion  libraries,  and  ineludes  the  ability  to  use  physies  models  for  energies  where 
tabular  data  are  not  available  {6).  MCNPX  treats  an  arbitrary  three-dimensional  eonfiguration  of 
materials  in  geometrie  eells.  Pointwise  eontinuous-energy  eross  seetion  data  for  seattering  in 
materials  are  used.  Souree  and  tally  struetures  with  extensive  statistieal  analysis  of  eonvergenee 
are  integral  to  the  output  parameters.  Rapid  eonvergenee  is  enabled  by  a  wide  variety  of 
variance  reduetion  methods.  MCNPX  ineludes  libraries  for  neutrons,  photons,  eleetrons,  protons 
and  photonuelear  interaetions.  “Mesh”  and  “radiography”  tallies  are  ineluded  for  2-  and 
3-dimensional  imaging  purposes.  An  auxiliary  program,  GRIDCONV,  converts  the  mesh  and 
radiography  tally  as  well  as  standard  metal-file  results  for  viewing  by  independent  graphies 
paekages. 

Monte  Carlo  algorithms  ealculate  a  statistieally  valid  result  of  the  most  likely  results  of  a  set  of 
operations  (7).  These  algorithms  are  useful  if  you  want  to  ask  what  will  be  the  outeome  of  a 
series  of  events.  This  is  partieularly  useful  in  eomplex  series  of  moleeular  motions  or  a 
sequential  series  of  events,  as  in  photon  seattering  through  a  material.  Eaeh  event  has  its  own 
probability  of  oeeurrenee.  It  ean  take  days  to  run,  in  order  to  test  as  many  events  and  outeomes 
as  possible.  These  algorithms  use  random  number  generators  to  ereate  the  spaee  for  variety  in 
the  eomplexity  and  follow  with  varianee  reduetion  teehniques  to  develop  statisties  that  give 
guidanee  as  to  the  value  or  weight  of  the  answer.  Credit  for  inventing  the  algorithm  often  goes 
to  Stanislaw  Ulam.  Unlike  Monte  Carlo  simulation,  deterministie  transport  methods  require  that 
the  phase-spaee  be  diseretized  in  spaee,  energy  and  angle,  introdueing  a  substantial  degree  of 
eomplieation  in  terms  of  algorithm  development,  problem  set-up,  and  requires  a  signifieant 
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amount  of  computer  memory  in  order  to  store  the  spatial,  angular  and  energy  dependent  flux 
information.  The  large  demand  placed  on  computer  memory  by  deterministic  transport  codes 
has  limited  their  use.  However,  with  the  rapid  increase  in  eomputer  processing  power,  memory 
capaeity,  and  reduction  in  unit  cost,  three-dimensional,  time-independent,  deterministic  transport 
codes  are  now  becoming  available  as  produetion  codes  (5). 

This  paper  describes  the  MCNPX  simulations  developed  to  determine  the  free  electron  density 
generated  in  bulk  SiC  layers.  The  P-partiele  souree  stimulation  modeled  in  the  numerieal 
simulation  is  the  speetrum  of  P  emitted  from  the  Sr90  isotope.  The  eleetron  density  gradient 
calculated  in  SiC  is  then  used  as  initial  conditions  for  a  drift-diffusion  numerical  model  of  a 
Sehottky  diode. 

Section  1  provides  the  motivation  for  our  efforts  and  discusses  some  of  the  applications  that  will 
benefit  from  the  effort.  Seetion  2  discusses  SiC  and  the  deviees  developed  for  this  application. 
Section  3  details  the  geometry  model  input  to  the  monte-carlo  code.  Sections  4  details  the 
significant  results  and  section  5  provides  conclusions  and  impact  of  this  part  of  the  larger  effort. 


2.  Materials 


Silieon  Carbide  (SiC)  is  a  robust  radiation-hard  solid-state  material  that  can  be  fabricated  into 
many  device  formats.  For  these  reasons  it  is  worth  pursuing  as  a  direct  energy  converter  (DEC) 
from  radiation  (alphas,  betas,  or  gammas)  into  free  eleetrons  that  can  be  used  as  a  source  of 
electrieal  current.  A  parametric  study  of  both  materials  and  deviee  geometry  is  underway. 
Materials  under  investigation  include  silicon  (Si  -2.3  g/cc),  SiC  (3.2  g/cc),  gallium  nitride  (GaN  - 
6.15  g/cc),  and  CVD  Diamond  (3.52  g/cc).  Device  geometries  are  not  identical  for  all  materials, 
however,  PIN  and  Sehottky  geometries  have  been  acquired.  Both  Si  and  SiC  Sehottky  diodes 
are  fabricated  at  the  Army  Researeh  Laboratory.  The  diodes  have  been  exposed  to  betas  emitted 
from  Strontium-90,  with  the  voltage  aeross  the  diode  measured  during  controlled  experiments. 
The  irradiation  of  these  deviees  from  strontium-90  will  be  the  focus  of  this  paper.  A  more 
complete  review  of  PIN  diode  structures  will  be  part  of  future  efforts. 

We  composed  the  results  of  numerieal  simulations  on  Al,  Si  and  SiC.  While  the  Si  and  SiC 
simulations  model  devices  were  built  and  measured  in  experiments,  the  purpose  of  the  aluminum 
(Al  -  2.7  g/ce)  is  1)  to  compare  to  the  Si  and  SiC,  and  2)  serve  as  a  reference  material  for 
electron  range  results.  Tables  of  eleetron  losses  have  been  tabulated  for  a  variety  of  materials 
including  Al  (P).  The  Berger  and  Seltzer  report  includes  tables  of  range  eonstants  (g/cm  )  as  a 
function  of  monochromatic  beta  energy.  The  incoming  beta  energy  used  as  a  eomparison  in  our 
calculations  is  150  keV  betas,  which  is  a  good  comparison  for  our  125  keV  (average)  beta 


2 


spectrum  used  in  the  MCNPX  calculations.  The  range  constant  supplied  in  section  3.2  (P)  must 
be  multiplied  by  the  density  of  the  material  used.  The  density  of  A1  is  shown  in  equation  1 . 

3.64x10'"^ - \ — ^  =  1.3xl0“"cm  =  130//w  (1) 

cm  Ug/cm 

The  range  of  the  aluminum  hit  by  the  constant  energy  150  keV  betas  is  shown  to  be  130  pm,  as 
shown  in  equation  1,  above,  which  is  a  good  comparison  to  the  150  pm  deep  model  used,  which 
will  be  described  in  more  detail  in  section  3.1. 


3.  Numerical  Simulations 


3.1  Geometry 

The  MCNPX  input  deck  that  defines  the  geometry  used  for  this  study  (see  appendix  A)  created 
thirty  rectangular  parallelepipeds  of  material  layers  in  a  row,  located  1  mm  from  a  0.5  mm  in 
diameter  beta  source  (see  figure  1).  The  physical  diode  package,  used  in  the  experiments, 
measures  5  mm  by  5  mm  by  0.5  mm.  However,  our  interest  was  to  see  the  number  of  electrons 
generated  in  the  top-most  layers  of  Schottky  diodes,  which  have  a  depletion  layer  of 
approximately  5  pm  thick.  Therefore,  the  layers  used  in  the  input  deck  measured  5  mm  by  5  mm 
in  the  plane  perpendicular  to  the  beta  source  and  5  pm  deep.  All  together,  the  thickness  of  the 
SiC  totaled  150  pm  to  ensure  a  complete  stopping  range  for  the  betas. 

It  is  essential  to  note  that  the  diode  interaction  with  the  electrical  circuit  cannot  be  simulated 
within  the  MCNPX  code.  Instead,  we  are  only  modeling  the  SiC  material  and  attempting  to 
calculate  how  many  free-electrons  are  being  generated  within  as  a  function  of  depth  in  the 
material.  A  separate  numerical  solid-state  code  will  be  used  in  an  additional  study  to  analyze  the 
electrical  response  of  the  diode  with  respect  to  an  external  circuit.  This  will  be  discussed  in 
subsequent  papers. 


3 


Figure  1.  Schematic  diagram  of  geometry  of  input  deck  of  MCNPX. 

3,2  Beta  Source 

The  beta  souree  coded  into  the  MCNPX  input  deck  (see  appendix  A)  is  an  attempt  to  model  a 
strontium-90  (Sr-90)  beta  emitter.  There  are  several  flaws  with  the  MCNPX  approximation, 
which  are  small,  but  necessary  to  recognize.  Strontium-90  emits  a  beta  spectrum  up  to  547  keV. 
However,  the  MCNPX  input  deck  was  simplified  to  an  input  spectrum  of  500  keV,  which  peaks 
at  125  keV  based  on  previous  calculations.  This  was  not  a  terrible  assumption  considering  that 
the  spectrum  of  electron  energies  emitted  from  the  isotope  is  usually  peaked  at  one  third  the 
maximum  energy  of  the  emitted  beta  {10).  The  beta  input  spectrum  used  to  simulate  the  Sr90 
isotope  is  shown  in  flgure  6.  Furthermore,  when  low  energy  electrons  enter  the  material,  fewer 
additional  electrons  are  generated,  though  the  energy  of  the  electrons  may  be  more  suitable  to  be 
captured  by  the  depletion  region  and  converted  into  usable  power  in  a  circuit. 


4.  Results 


Contour  plots  of  electron  number  density  (#/cc),  shown  in  figure  2,  were  generated  for  each  of 
the  three  materials  being  evaluated;  SiC,  Al,  and  Si.  The  calculation  using  MCNPX  gives 
results  of  energy  deposited  with  units  of  MeV/g.  The  number  density  is  calculated  from  these 
original  results  by  a)  multiplying  by  material  density  and  b)  dividing  by  the  average  particle 
energy  (125  keV).  The  number  density  data  was  calcualted  by  recording  MCNPX  mesh 
diagnostic  over  a  volume  that  spanned  5  mm  by  0.15  mm  by  0.15  mm.  The  SiC  layer  volume 
modeled  in  the  MCNPX  input  deck  occupies  a  total  of  5  mm  by  5  mm  by  0.15  mm.  In  order  to 
get  a  number  density  at  the  center  of  the  beam,  the  mesh  diagnostic  was  implemented  to 


4 


calculate  the  energy  deposited  over  a  much  smaller  volume  (at  the  beam  center)  than  the  entire 
layer.  The  binning  of  the  mesh  data  cell  was  defined  to  be  250  pm  by  8  pm  by  150  pm.  By 
taking  a  smaller  mesh  size,  a  more  accurate  description  was  obtained.  The  results  of  the  contour 
plots  show  that  the  number  density  at  the  center  of  the  beam  is  27  k  electrons/cc/incident 
electron. 


Figure  2.  Number  density  contour  plot  for  (from  left  to  right)  (a)  SiC,  (b)  Al,  and 
(c)  Si.  The  view  is  5  mm  by  0.15  mm,  and  the  slice  is  0.15  mm  thick, 
representing  30-5  pm  partial-layers.  The  binning  is  stored  in  such  a  way 
that  each  data  point  measures  250  pm  by  8  pm  by  150  pm.  Closer 
inspection  shows  that  the  peak  number  density  decreases  by  approximately 
20%  in  each  case. 

Number  density  plots  (figure  3)  show  a  1 -dimensional  horizontal  out  aoross  the  oontour  plots  of 
figure  2.  They  ourves  show  more  quantitatively,  the  number  density  values  along  the  oenter-line 
(beam-line)  in  1 -dimension,  oompared  to  what  is  piotured  in  2-D  in  the  oontour  plots  of  figure  2. 
With  the  addition  of  a  trend  line  to  eaoh  data  set,  the  number  density  is  shown  to  deorease 
exponentially  as  eleotrons  propagate  the  layers  of  SiC.  Variation  from  the  line  may  result  from 
poor  statistios,  whioh,  though  satisfaotory,  oould  be  improved  with  more  histories  numerioally 
oaloulated  (and  longer  run  times).  The  statistical  result  of  the  mesh  oaloulation  is  well  below 
10%  which  provides  a  high  degree  of  eertainty. 
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(a) 


(b) 


(c) 


Figure  3.  Number  density  plot  for  (from  top  to  bottom)  (a)  SiC,  (b)  Al,  and 
(c)  Si.  The  data  is  taken  along  the  eentral  line  of  the  layers,  where 
eaeh  point  represents  250  pm  by  8  pm  by  150  pm.  The  data  is  shown 
with  exponential  fits.  The  point  shown  in  the  lower  right  of  eaeh 
plot  is  the  exponential  fit  to  the  order  of  magnitude  fall  off  from  the 
starting  point. 
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We  expected  that  the  number  of  free  electrons  generated  in  the  materials  would  be  correlated  to 
material  density.  They  are  not,  however,  linearly  dependent  on  material  density.  With  only 
three  data  points  to  compare,  it  is  left  for  a  future  study,  to  better  understand  this  result.  One  way 
to  understand  this  discrepancy,  is  that  not  only  was  the  material  density  changing,  but  so  were 
the  atomic  scatterers.  The  density  varied  during  the  three  calculations  from  2.3,  2.7,  and 
3.2  g/cc,  while  the  atomic  number  of  the  materials  also  changed  from  Si(atomic  mass  28.1), 
Al(atomic  mass  26.9),  and  SiC.  A  linear  approximation  is  still  attempted  to  fit  the  density  of  free 
electrons  generated  by  the  Sr90  betas  with  changing  density  of  the  DEC  material.  This  linear  fit 
to  the  data  is  shown  below  in  figure  4. 


The  electron  range  of  each  of  the  three  DEC  materials  is  plotted  (figure  4)  as  a  function  of  the 
DEC  material  density.  The  electron  range  here  is  defined  as  10%  of  initial  number  density  at  the 
surface.  Before  extrapolation,  all  the  materials  looked  like  they  would  have  a  similar  fall  off 
point,  though  closer  inspection  revealed  that  the  number  density  has  a  linear  relationship  with 
respect  to  material  density,  just  as  the  peak  number  density  does.  However,  more  data  would  be 
needed  for  a  better  evaluation. 
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Figure  4.  Peak  eleetron  number  density  (a)  data  gathered  from  the  Number  Density  Plots.  The 

exponential  fit  was  was  extrapolated  to  zero,  and  the  zero  points,  or  starting  points  of  the 

materials  were  plotted  with  respeet  to  the  density  of  the  material.  Electron  range 

(b)  data  shows  the  point  at  which  the  density  is  10%  of  the  value  at  the  surface  of  the  material. 


The  spectrum  of  the  electrons  depositing  energy  in  a  layer  is  shown  in  figure  5.  This  graph 
clearly  shows  that  the  lower  energy  electrons  are  responsible  for  depositing  most  of  the  energy  in 
the  system.  The  spectrum  shown  is  unlike  the  shape  of  the  input  Sr90  spectrum.  The  beta 
spectrum  bombarding  the  SiC  layers  is  shown  in  figure  6.  The  comparison  of  input  beta 
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spectrum,  and  spectrum  of  electrons  depositing  energy  in  the  material  shows  that  the  higher 
energy  betas  tend  to  travel  through  the  material  without  leaving  much  energy  behind.  This 
would  lead  us  to  believe  that  the  scattered,  lower-energy  (<100  keV)  electrons  are  more 
significant  in  providing  power  to  an  external  circuit. 


Figure  5.  Energy  speetrum  of  number  density  of  eleetrons  depositing  energy  into  layers  1  through  3  of  the 
5  pm  thiek  SiC  layers. 
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Figure  6.  Input  beta  spectrum  from  Strontium  90  is  peaked  at  125  keV.  This  spectrum  is  the  input  modeled 
in  the  MCNPX  simulation. 


5.  Conclusions 


The  number  density  result  calculated  in  MCNPX  is  the  number  of  electrons  generated  for  each 
incident  electron.  The  Sr90  source  used  in  our  experiments  is  physically  located  at  NSWC- 
Carderock.  It  has  an  activity  of  7mCi.  This  corresponds  to  25e7  beta  emissions  per  second  (Bq). 
The  resulting  number  density  of  electrons  generated  as  a  result  of  our  experimental  setup  is  then 
calculated  in  the  SiC  converter  experiments  to  be  675el0/cc/sec. 

The  range  for  betas  emitted  from  Strontium  90  has  been  modeled  using  MCNPX.  The  results 
show  that  150  pm  of  SiC  stops  all  electrons  emitted.  A  parametric  study  of  target  material  was 
performed  using  MCNPX.  Input  decks  for  this  numerical  investigation  differ  only  by  material 
(Al,  Si,  and  SiC)  and  their  corresponding  naturally  occurring  density.  These  results  compare 
well  with  values  generated  by  Berger  and  Setter  (5)  in  which  they  tabulated  the  values  of 
stopping  range  in  Al,  Si  and  Cu. 

Diagnostics  available  in  the  MCNPX  simulations  were  compared  for  consistency.  Surface  flux 
results  (tally  2)  compare  well  with  volume  flux  results  (tally  4),  in  the  limit  as  the  thickness  of 
the  volume  decreased  to  a  surface. 
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When  calculating  the  energy  deposited  in  the  material,  tally  6  generates  a  result  over  a  complete 
cell.  When  a  more  localized  or  discretized  value  is  required  (i.e.,  the  value  at  the  center  of  the 
beam),  using  the  mesh  command  and  identifying  smaller  effective  volume  regions  works  well 
and  compares  accurately  to  the  averaged  value  over  the  larger  volume  of  the  cell.  The  value 
calculated  for  energy  deposited  per  incident  electron  per  gram  of  material  is  22.6  MeV/g.  The 
layer  dimension  is  5  mm  x  5  mm  x  5  pm,  or  1.25e-4  cc.  The  density  of  SIC  is  3.2  g/cc 
corresponding  to  a  layer  mass  of  4e-4  g.  The  resulting  energy  deposited  is  9.04  keV.  This 
represents  7%  of  the  peak  incident-electron-energy  (125  keV),  deposited  in  a  5  um  layer  of  SIC. 

In  a  prototype  isotope  battery  design,  foils  will  be  sandwiched  in-between  two  SiC  DEC  devices. 
In  order  to  better  understand  the  efficiency  that  would  be  expected  from  an  isotope  foil  above  the 
SiC  DEC  layer,  we  performed  a  broad  beam  simulation.  The  incident  beam  radius  was  varied 
[0.5,  1,  2,  5  mm].  Beams  diameters,  less  than  2  mm,  generated  electron  showers  completely 
within  the  5  mm  side-length  of  the  devices  under  test.  The  energy  deposited  in  the  first  5  pm 
layer  was  calculated  to  be  22.6  MeV/g.  2.26%  of  the  incident  energy  is  deposited  in  the  first 
layer.  This  result  is  identical  to  the  previous  pencil  beam  simulations. 

These  nuclear  scattering  calculations  do  not  take  into  account  the  electrical  circuit  properties  of  a 
diode  (i.e.,  source  impedance  feeding  a  load  impedance).  The  numerical  simulations  described 
above  are  meant  to  determine  an  initial  electron  density,  or  density  of  carriers  for  initial 
conditions  into  a  circuit-based  model.  The  initial  conditions  are  used  in  a  study,  already 
underway,  to  model  the  drift  and  diffusion  in  diode  structures.  The  results  will  be  reported  in  our 
next  technical  report  of  the  investigation. 

The  follow-on  circuit  modeling  is  expected  to  provide  current  and  voltage  numerical  results  and 
insight  into  the  electron  flow  within  the  semiconductor  materials  that  convert  the  emitted-beta 
energy  into  electrical  current  for  a  long-lived  power  source.  The  results  of  the  circuit  models 
being  developed  will  be  compared  to  measurements.  With  electron  density  initial  conditions  (the 
result  of  this  investigation),  a  circuit  model,  and  measurements  in-hand,  an  optimized  DEC  for  a 
nuclear  battery  design  is  expected  as  a  result  in  the  near  future. 
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Appendix  A.  MCNPX  Input  Deck 


The  input  deck  shown  below  is  one  of  the  series  of  numerical  simulations  evaluated. 


sicm3 

Model  of  Silicon  Carbide  with  Sr  source 
c  (MCNPX  scaled  in  cm) 
c 


c  This  is  a  model  of  the  way  electrons  would  be  generated 
c  in  the  SiC  after  being  hit  by  electrons  from  a  strontium  source. 

c 

c  + . 

c| 

. + 

c| 

c| 

C  + . 

Cell  Cards 

1 

. + 

c  #  mat  density  surface  data 

cSiC 

1  3  -3.2 

-1 

imp:e  1  u=l 

$SiC  layer  closest  to  source 

2  3  -3.2 

-2 

imp:e  1  u=l 

$2nd  SiC  layer 

3  3  -3.2 

-3 

imp:e  1  u=l 

$3rd  SiC  layer 

4  3  -3.2 

-4 

imp:e  1  u=l 

$4th  SiC  layer 

5  3  -3.2 

-5 

imp:e  1  u=l 

$5th  SiC  layer 

6  3  -3.2 

-6 

imp:e  1  u=l 

$6th  SiC  layer 

7  3  -3.2 

-7 

imp:e  1  u=l 

$7th  SiC  layer 

8  3  -3.2 

-8 

imp:e  1  u=l 

$8th  SiC  layer 

9  3  -3.2 

-9 

imp:e  1  u=l 

$9th  SiC  layer 

10  3  -3.2 

-10 

imp:e  1  u=l 

$10*  SiC  layer 

11  3  -3.2 

-11 

imp:e  1  u=l 

$1 1th  SiC  layer 

12  3  -3.2 

-12 

imp:e  1  u=l 

$12th  SiC  layer 

13  3  -3.2 

-13 

imp:e  1  u=l 

$13th  SiC  layer 

14  3  -3.2 

-14 

imp:e  1  u=l 

$14th  SiC  layer 

15  3  -3.2 

-15 

imp:e  1  u=l 

$15th  SiC  layer 

16  3  -3.2 

-16 

imp:e  1  u=l 

$16th  SiC  layer 

17  3  -3.2 

-17 

imp:e  1  u=l 

$17th  SiC  layer 

18  3  -3.2 

-18 

imp:e  1  u=l 

$18*  SiC  layer 

19  3  -3.2 

-19 

imp:e  1  u=l 

$19*  SiC  layer 

20  3  -3.2 

-20 

imp:e  1  u=l 

$20*  SiC  layer 

21  3  -3.2 

-21 

imp:e  1  u=l 

$21th  SiC  layer 

22  3  -3.2 

-22 

imp:e  1  u=l 

$22th  SiC  layer 

23  3  -3.2 

-23 

imp:e  1  u=l 

$23th  SiC  layer 

24  3  -3.2 

-24 

imp:e  1  u=l 

$24th  SiC  layer 

25  3  -3.2 

-25 

imp:e  1  u=l 

$25th  SiC  layer 

26  3  -3.2 

-26 

imp:e  1  u=l 

$26th  SiC  layer 

27  3  -3.2 

-27 

imp:e  1  u=l 

$27th  SiC  layer 

28  3  -3.2 

-28 

imp:e  1  u=l 

$28*  SiC  layer 

29  3  -3.2 

-29 

imp:e  1  u=l 

$29*  SiC  layer 

30  3  -3.2 

c 

-30 

imp:e  1  u=l 

$30*  SiC  layer 

c  air  encasing 

31  2  -0.00129  -31  1  2345  6  7  8  9  10  11  12  13  14  15  16  17  18  19  & 
20  21  22  23  24  25  26  27  28  29  30  imp:e  1  u-1 

c 

32  0  1  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16  17  & 

18  19  20  21  22  23  24  25  26  27  28  29  30  imp:e  1  u-1 


33  0  31  imp:e  1  u=l 

c  room 

34  2  -0.00129  -32  imp:e  1  fill-1 

35  0  32  imp:e0  Soutside  world 

c  + . + 

c| 

c  I  Surface  Cards 

c| 

c  + . + 


c  #  type  params 
cSiC 

1  rpp  708.75  709.25  493.9995  494  117.25  117.75  $SiC  layer  closest  to  source 

2  rpp  708.75  709.25  493.9990  493.9995  117.25  117.75  $2tid  SIC  layer 

3  rpp  708.75  709.25  493.9985  493.9990  117.25  117.75  $3rd  SiC  layer 

4  rpp  708.75  709.25  493.9980  493.9985  117.25  117.75  $4th  SiC  layer 

5  rpp  708.75  709.25  493.9975  493.9980  117.25  117.75  $5th  SiC  layer 


6  rpp  708.75  709.25  493.9970  493.9975  117.25  117.75  $6th  SiC  layer 

7  rpp  708.75  709.25  493.9965  493.9970  117.25  117.75  $7th  SiC  layer 

8  rpp  708.75  709.25  493.9960  493.9965  117.25  117.75  $8th  SiC  layer 

9  rpp  708.75  709.25  493.9955  493.9960  117.25  117.75  $9th  SiC  layer 

10  rpp  708.75  709.25  493.9950  493.9955  117.25  117.75  SlOth  SiC  layer 

11  rpp  708.75  709.25  493.9945  493.9950  117.25  117.75  $1  Ith  SiC  layer 

12  rpp  708.75  709.25  493.9940  493.9945  117.25  117.75  $12*  SiC  layer 

13  rpp  708.75  709.25  493.9935  493.9940  117.25  117.75  $13*  SiC  layer 

14  rpp  708.75  709.25  493.9930  493.9935  117.25  117.75  $14tliSiClayer 

15  rpp  708.75  709.25  493.9925  493.9930  117.25  117.75  $15tliSiClayer 

16  rpp  708.75  709.25  493.9920  493.9925  117.25  117.75  $16tliSiClayer 

17  rpp  708.75  709.25  493.9915  493.9920  117.25  117.75  $17tliSiClayer 

18  rpp  708.75  709.25  493.9910  493.9915  117.25  117.75  $18tliSiClayer 

19  rpp  708.75  709.25  493.9905  493.9910  117.25  117.75  $19*SiClayer 

20  rpp  708.75  709.25  493.9900  493.9905  117.25  117.75  $20*  SiC  layer 

21  rpp  708.75  709.25  493.9895  493.9900  117.25  117.75  $21*  SiC  layer 

22  rpp  708.75  709.25  493.9890  493.9895  117.25  117.75  $22*  SiC  layer 

23  rpp  708.75  709.25  493.9885  493.9890  117.25  117.75  $23*  SiC  layer 

24  rpp  708.75  709.25  493.9880  493.9885  117.25  117.75  $24*  SiC  layer 

25  rpp  708.75  709.25  493.9875  493.9880  117.25  117.75  $25*  SiC  layer 

26  rpp  708.75  709.25  493.9870  493.9875  117.25  117.75  $26*  SiC  layer 

27  rpp  708.75  709.25  493.9865  493.9870  117.25  117.75  $27*  SiC  layer 

28  rpp  708.75  709.25  493.9860  493.9865  117.25  117.75  $28*  SiC  layer 

29  rpp  708.75  709.25  493.9855  493.9860  117.25  117.75  $29*  SiC  layer 

30  rpp  708.75  709.25  493.9850  493.9855  117.25  117.75  $30*  SiC  layer 

c  lead  block  to  encase  silicon  carbide 

31  rpp  708.5  709.5  493.850  494  117  118 

c  room  dimensions 

32  rpp  0  1419  0  1984  0  466 

c  + . + 

c| 

c  I  Material  Cards  | 

c|  I 

c  + . + 

c  m#  isotope  percent  & 

c  ml  82000  1.0  Slead  11.34  g/cc 

m2  8016  .3  & 

7014  .7  $air  0.00129  g/cc 

m3  6012  .43  & 

14028  .57  $SiC  3.2  g/cc 

c  m4  14028  1.0  $Si  2.3  g/cc 

c  ml3  6012  1.0  $C  2.3  g/cc 

c  m5  90000  1.0  $Sr 

c 

c  + . + 

c| 

c  I  Source  Definition  | 

c|  I 

c  + . + 

sdef  & 
dir^l  & 
erg^  d2  & 
vec^O  -1  0  & 
par=e  & 
rad^dl  & 
pos-709  497  117.5 
c 

sil  0  .05 
spl  -21 
c 

c  energy  distribution  taken  from  now  strontium  (old:  280kv-8ma-12may-aft.spe) 
si2  0.000916  0.00549  0.0412  0.0833  0.125  0.167  0.208  0.25  0.291  0.333  & 
0.375  0.417  0.459  0.5 

sp2  0  1000  3405  4692  4768  3900  2775  2234  1725  1331  & 

1007  764  556  378 
c 

c  si3  h  .707  1 
c  sp3  d  0  1 
c 
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c  + . + 

C|  I 

c  I  Data  Cards  | 

c|  I 

c  + . + 

mode  e 

nps  10000  $10Mh-9hrs 
c 

c  + . + 

c| 

c  I  Tallies  | 

c|  I 

c  + . + 

c 

c 

fc2  test 
f2:e  1.3 


e2  0  98i  .51 
c 

c  fc4  electron  flux  averaged  over  the  cell  in  particles/cm'^2 
c  f4:e  1  2  3 
ce4  0  98i  .51 


cora21  708.75  18i  709.25 
corb21  493.9995  18i  494 
corc21  117.25  117.75 

c 

rmesh3 1  :e  flux  pedep  popul  dose 
c  layer  1  SiC 

cora31  708.75  709.25 

corb31  493.9995  18i  494 
corc31  117.25  18i  117.75 
c 
c 

rmesh41:e  flux  pedep  popul  dose 
c  layer  2  SiC 

cora41  708.75  18i  709.25 
corb41  493.999  493.9995 

corc41  117.25  18i  117.75 
c 

rmesh5 1  :e  flux  pedep  popul  dose 
c  layer  2  SiC 

coraSl  708.75  18i  709.25 
corb51  493.999  18i  493.9995 
corc51  117.25  117.75 


c 

c 

c  this  tally  may  not  work  {or  be  correct)  because  it  is  with  electrons 
c  fc6  electron  energy  deposition  averaged  over  the  cell  in  MeV/g 
c  f6:e  1  2  3 
ce6  0  98i  .51 


c 

rmesh61:e  flux  pedep  popul  dose 
c  layer  2  SiC 

cora61  708.75  709.25 

corb61  493.999  18i  493.9995 
corc61  117.25  18i  117.75 


c 

c  fc8  electron  energy  deposition  in  MeV 
c*fS:e  123 
ce8  0  le-5  97i  .51 


c 

c  + . + 

c| 

c  I  Mesh  I 

c|  I 

c  + . + 

c 


tmesh 

c 

c  (a)radial  meshes  (b)zaxis(symmetry)  (c)angles  counterclockwise 
c 

c  flux=#/cni2  $dose=rem/hr  $pedep=MeV/cm3 
c 

c  flux:  flux 

c  pedep:  total  energy  deposition/unit  vol 
c  dose:  dose  over  cross  planes 


c 

c 

rmesh71:e  flux  pedep  popul  dose 
c  layer  3  SiC 

cora71  708.75  18i  709.25 
corb71  493.9985  493.999 

corc71  117.25  18i  117.75 
c 

rmesh81:e  flux  pedep  popul  dose 
c  layer  3  SiC 

cora81  708.75  18i  709.25 
corbSl  493.9985  18i  493.999 
corcSl  117.25  117.75 

c 

rmesh91:e  flux  pedep  popul  dose 
c  layer  3  SiC 

cora91  708.75  709.25 

corb91  493.9985  18i  493.999 
corc91  117.25  18i  117.75 
c 


c 

rmeshl  :e  flux  pedep  popul  dose 
c layer  1  SiC 

coral  708.75  18i  709.25 
corbl  493.9995  494 

cord  117.25  18i  117.75 
c 

rmesh21:e  flux  pedep  popul  dose 
c  layer  1  SiC 


c 

rmeshl01:e  flux  pedep  popul  dose 
c  all  SiC  layers  (30) 
coralOl  708.75  18i  709.25 
corblOl  493.985  25i  494 
corclOl  117.4925  117.5075 

c 

endmd 
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